基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

导流槽构型对火箭起飞时流场和声场的影响

为了研究四喷管运载火箭起飞时火箭周围噪声环境问题,建立燃气/空气双组分的可压缩流动模型,采用2阶Roe格式、SAS(scale-adaptive simulation)湍流模型和声学类比积分法Ffowcs-Williams Hawkings(FW-H)求解三维Navier-Stokes方程。以单机火箭的噪声问题为对象开展数值模拟,并将噪声数值计算结果与试验数据对比,误差在3dB以内(相对误差小于1.6%),验证预测噪声方法的有效性,进而研究四喷管运载火箭起飞阶段采用不同导流槽构型对箭体舱段区域噪声环境的影响,结果表明:在相同噪声接收点处,单侧导流槽对应的总声压级比双侧导流槽大,两者的总声压级(OASPL)之差最大为10.7dB。另外,当采用单侧导流槽时,沿周向接收点的噪声总声压相对单侧导流槽中心截面呈对称分布,而且沿导流出口方向逐渐增大。所建立的噪声数值方法为大推力捆绑运载火箭舱段的噪声环境预测及其控制提供一定参考。     

大涵道比涡扇发动机TPS短舱低速气动特性分析

为了评估民机低速带动力试验时进排气效应的影响,选取大涵道比发动机涡轮动力模拟器(TPS)短舱和真实发动机短舱作为研究对象,采用数值模拟方法对其起飞、进近状态的低速气动特性进行对比分析。结果表明:由于TPS流量低于真实发动机需求,其唇口、外罩流场特征和真实发动机短舱有所不同,阻力特性也有差别;在进气道处于亚临界状态时,TPS短舱阻力系数比真实短舱大了约1.7个阻力单位,又由于唇口当地气流攻角更大,使得TPS短舱失速攻角相对降低了约1.0°;当进气道工作于超临界状态时,TPS短舱虽然也可以反映真实短舱的流动特性,但由于捕获流管收缩情况和气流驻点随攻角的变化,使得在0°~20°攻角时TPS短舱的阻力系数高于真实短舱,而在20°~30°攻角时其阻力系数略低,差量最大约为1.8个阻力单位。对于研究的大涵道比发动机,未经唇口及外罩修正的TPS短舱其低速气动特性基本可以反映真实进排气效应的影响,但在气动特性分析中可以考虑进一步修正进气效应的影响。     

乙醇-空气稀液雾值班火焰的大涡模拟研究

为加深对气液两相湍流燃烧现象的认知,检测火焰面模型在液雾燃烧中的适用性,本文在欧拉-拉格朗日架构下使用火焰面/反应进度变量模型(FPV)数值模拟了乙醇-空气稀液雾值班火焰。欧拉坐标系下的气相湍流场使用大涡模拟方法模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑了相间质量、动量、能量交换。燃烧模型中采用碳元素定义混合物分数,在混合物分数方程源项中体现液相对燃烧模型的影响。模拟得到的气相温度分布和液相统计值均和实验数据较好吻合,验证了该燃烧模型对稀液雾扩散类型火焰的适用性。分析瞬时图发现,该稀液雾火焰的最高燃烧温度往往位于当量混合物分数附近,在出口下游20倍直径处火焰完全点着,此处上游FPV模型能给出局部点火熄火现象。蒸发作用在剪切层和点火区域较强,而液滴和火焰作用较弱,单一液滴很少被火焰包围。     

跨声速压气机多圆柱孔式处理机匣设计与扩稳机理研究

基于对跨声速轴流压气机内部失稳触发机制的认识,设计了一种新型多圆柱孔式处理机匣结构。对带该处理机匣的跨声速轴流压气机转子Rotor-37的内部流动进行了全三维非定常数值模拟,结果表明该新型多圆柱孔式处理机匣结构可使压气机的综合裕度提高6.5%,而最高效率点效率仅降低0.19%。对多圆柱孔式处理机匣的扩稳机理进行了详细分析,结果表明处理机匣能有效消除低速流体团引起的通道阻塞现象,并将低速阻塞流体团抽吸进入处理机匣容腔,然后从前端喷口预旋喷出以抑制叶顶区域的流动分离,从而有效提高跨声速轴流压气机的失速裕度。     

基于环控物理仿真平台的余度模型设计

在环控物理仿真平台上,基于MATLAB和Vxworks环境,建立了飞机环控系统余度模型。详细论述了环控系统故障模型、机电管理计算机模型。此外,对模型的有效性、可靠性进行了验证。通过试验验证,表明模型结构正确、能够满足环控系统余度功能要求。     

基于模型的系统工程在作动器领域的应用研究

以飞机作动器产品为研究对象,通过条目化需求分析,将产品设计过程进行逐级分解。利用AMESim、AN-SYS系列软件等仿真工具,建立从部件到产品的各级仿真模型,并开展仿真分析工作,通过对仿真分析过程和结果的分析,探索了基于模型的系统工程方法的工作流程,总结了飞机液压作动器基于模型的系统工程方法的应用基础。     

7150铝合金热变形行为及微观组织

采用Gleeble热模拟机进行热压缩实验,研究7150铝合金在变形温度为300~450℃、应变速率为0.01~10s-1条件下的变形行为,采用Zener-Hollomon参数法构建合金高温塑性变形本构方程,并对变形后的微观组织进行分析。研究表明:7150铝合金的流变应力随应变速率增大而增大,随变形温度增大而降低。该合金热压缩变形的流变应力行为可用双曲正弦形式的本构方程描述,其参数A为4.161×1014s-1,α为0.01956 MPa-1,n为5.14336,热变形激活能Q为229.7531k J/mol。随着温度升高和应变速率降低,动态再结晶逐渐取代动态回复成为合金的主要软化机制。     

多轮多支柱起落架飞机操纵前轮转弯特性分析

在考虑机体弹性的前提下,采用Fiala轮胎理论,建立了多轮多支柱式起落架飞机全机动力学模型,进行飞机操纵前轮转弯分析。通过仿真,给出了典型前轮操纵角下,各轮胎上的径向及侧向载荷。并以轮胎出现侧滑为临界条件,研究了不同前轮操纵角下允许的飞机最大滑行速度,给出飞机操纵前轮转弯包线。最后研究了飞机180°转弯时重心及轮胎的运动轨迹,得出在保证操纵安全的前提下,飞机可在50 m宽的跑道上实现180°转弯。